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Abstract 

We study scalar li-variate subdivision schemes, with dilation matrix 21, sat- 
isfying the sum rules of order k. Using the results of MoUer and Saner, stated 
for general expanding dilation matrices, we characterize the structure of the 
mask symbols of such schemes by showing that they must be linear combina- 
tions of shifted box spline generators of some quotient polynomial ideal. The 
directions of the corresponding box splines are columns of certain unimodular 
matrices. The quotient ideal is determined by the given order of the sum rules 
or, equivalently, by the order of the zero conditions. 

The results presented in this paper open a way to a systematic study of 
subdivision schemes, since box spline subdivisions turn out to be the building 
blocks of any reasonable multivariate subdivision scheme. 

As in the univariate case, the characterization we give is the proper way 
of matching the smoothness of the box spline building blocks with the order of 
polynomial reproduction of the corresponding subdivision scheme. However, due 
to the interaction of the building blocks, convergence and smoothness properties 
may change, if several convergent schemes are combined. 

The results are illustrated with several examples. 
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Introduction 

Subdivision schemes are efficient iterative procedures for generating finer 
and finer grids of points in M'^ and are used to design smooth curves or surfaces. 
Starting with some initial grid of points, a binary scalar subdivision scheme 
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computes the coordinates of the finer grid points d'^'""'"^^, inheriting the topology 
of the coarser ones, via local averages 

d(.+i)=5^dW^ J2 d^'-)a„_2/3, r>0. 

These averaging rules depend on the coefficients of the corresponding subdivision 
mask 

) cjd which we assume to be some finitely supported sequence of 
real numbers. If the mask is chosen appropriately, the grids with the vertices 
d^'^\ as r goes to infinity, can be interpreted as approximations of the values 
of a limiting curve or surface. The locality of the method and its algorithmic 
simplicity ensure that the subdivision recursion is fast, efficient, and easy to 
implement. These features and the connection between subdivision and other 
multiresolution methods have led to an increasing popularity of subdivision in 
computer graphics, in computer aided geometrical design, and in wavelet and 
frame constructions. For more details on subdivision we refer the interested 
reader to the pioneering work in 0, or to the more recent survey [l9[ and 
the references therein. For basic details on wavelet and frame constructions see, 
e.g., EE, Hi- 



In general, the limit of a convergent subdivision scheme is not known ana- 
lytically. Nevertheless, various analytical properties of the limit can be read off 
the mask symbol 



ci{z) = J2 ^ ^ (C\{0}) 



d 



which — due to the finite support of the mask — is a d-variate Laurent polynomial. 
In the univariate setting, the properties of convergent subdivision schemes are 
well-understood. The symbol a,{z) of any convergent univariate subdivision 
scheme possesses the Factorization Property: it can be written as 

a(^) = 2.a(z). with a{l) = l, 

for some integer fc > 1, with the Laurent polynomial (y{z) satisfying some addi- 
tional properties to ensure the convergence and smoothness of the subdivision 
limit, see [l^. [ioj. Since the coefficients of a sum to 1, this Factorization Prop- 
erty tells us that convergence of the subdivision scheme implies that its mask 
can be written as an affine combination of shifted versions of B-spline symbols 
Bk{z) = (1 -I- z)'^/2'^~^. Thus, any reasonable univariate subdivision scheme 
uses B-spline symbols as its building blocks, a fact which is ver^crucial for 
designing new efficient subdivision processes, see e.g., [HE El, 



In this paper we are going to present results which aim at replacing the 
Factorization Property of univariate schemes by a Decomposition Property of 
multivariate schemes, see Theorems A-C below. Theorem C deals with the case 
d = 2 and is crucial for studying the properties of subdivision schemes. It tells 
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us that the mask symbol of any reasonable bivariate subdivision scheme can be 
decomposed as 

a,{zi,Z2) = 4:- ^ >^a,p.,-f ■ Cra.^l3,j{zi,Z2) ■ B*jj^^{zi,Z2) , 



where ^a.i3,j = 1, and the symbols cTq,^^ Z2) are Laurent polynomials 
normalized by 0"^^^^^ (1,1) = 1. Here, 



B* 



are the normalized mask symbols of three-directional box splines, fc refers to the 
order of polynomial reproduction of the subdivision operator — in the wavelet 
literature also denoted as order of accuracy of the mask — and the list I/c is as 
in Theorem 13.21 The advantage of such a decomposition result is obvious, not 
only for the purpose of classifying the zoo of subdivision schemes, but also as a 
starting point for constructing new schemes or for enhancing the properties of 
the schemes according to specific requirements. 

We believe that the use of the building blocks ^ in the above decompo- 
sition is the appropriate generalization to the bivariate case of the normalized 
univariate B-spline symbols, since B^ ^ ^ correspond to a class of well-known 
bivariate spline functions whose order of global smoothness and whose order of 
polynomial reproduction match in the same way as in the univariate case, see 
the notes at the end of Section [3Jl We consider, therefore, this paper to be also 
an interesting contribution to the box spline literature. 

The methods employed in this paper are of algebraic nature, as well as other 
results on subdivision dealing with polynomial reproduction, see 0, Chapter 6]. 
It is well-known that the necessary conditions for the convergence of a subdivi- 
sion scheme are equivalent to the fact that the mask symbol a(2) is a properly 
normalized element of a certain ring of Laurent polynomials: the symbol must 
belong to the ideal T, or to its power , depending on the order k of the poly- 
nomial reproduction of the scheme. For definitions of I*^, fc G N, see (|2.1I) and 
(|1.6p . or see the zero conditions of order fc in (|1.7p . This alg ebraic property of 
&{z) has been studied in detail in the papers by Sauer 3^, 3lf and by Moller and 



Saner [28|, which also motivated us to consider these ideals once again. Note 
that the Factorization Property says that &.(z) belongs to the principal ideal 
generated by some B-spline symbol, while the Decomposition Property tells us 
that a(2:) is in the ideal generated by a particular set of box spline symbols. 
Moller and Sauer have given other sets of generators for these ideals, putting 
the emphasis on the algebraic properties of these generators. We consider it one 
of our main achievements that we were able to relate the generators of the ide- 
als to a well-known class of spline functions, see [31. Box spline subdivision 
has been studied thoroughly in the literature, see where such box spline 

schemes also appear under a different name. In the three-directional case, for 
example, their mask symbols are given by 4 • i?^^ ^(-^ij •2:2). 



3 



Starting with Section 2, we restrict ourselves to working with the ideals of 
the ring of d-variate polynomials instead of the ring of Laurent polynomials. 
Note that any Laurent polynomial can be shifted to produce a polynomial. It is 
done by multiplying the Laurent polynomial with a factor 2;°, i.e., a unit in the 
ring of Laurent polynomials. This results in a shift of the support of the mask. 
Such a shift does neither affect the convergence and regularity of the scheme, 
nor does it change the zero conditions of any order. We, therefore, assume that 
the mask is supported in Nq. 

The paper is organized as follows: in Section 1 we introduce some nota- 
tion and background on scalar subdivision schemes with dilation matrix 21. In 
Section 2 we show how the zero conditions of order fc = 1 on the mask symbol 
&{z) determine the structure of the symbol of any convergent binary subdivision 
scheme. This result is given in Theorem A and is a crucial step toward under- 
standing the case fc > 1. As in the multivariate case the Factorization Property 
is replaced by the Decomposition Property, it is important to reduce the number 
of the terms of this decomposition or, equivalently, to minimize the number of 
the generators for I, see Theorem 12. 6[ Theorem 12.111 Section 3 contains our 
main results. Theorems B and C, dealing with powers of I. The fact that the 
ideal I*^ is generated by the appropriate products of the elements of I makes 
the proof of Theorem B straightforward, if one is not interested in reducing the 
number of terms of the corresponding decomposition. We address the latter 
issue only in the case d = 2, as this case is of special interest in subdivision. 
The corresponding result is stated in Theorem C. In Section 4 we illustrate the 
result of Theorem C with several bivariate examples. 

This paper is an extended version of the technical report Q which can be 
obtained from the first author. 

1. Background on subdivision and notation 

A scalar d-variate subdivision scheme is given by a scalar Z'^-indexed se- 
quence a = (aa)aez<'7 the so-called mask, defining the subdivision operator <Sa 
on data sequences d = [Aot)aei.'^ ^ £(Z'*) as follows: 



We assume that the mask is finite, i.e., only finitely many coefficients 

t^ot are 

non-zero. 

In our study we use the following symbol notation. For a finitely supported 
sequence c — (Cq,)^^^!!, its symbol is given by the Laurent polynomial 




(1.1) 



c{z) = ^ c„ 2;' 



with z = {zi, . . . , Zd) e (C \ {0}) and, in the multi-index notation, 
z"" = zj'iz^^ - ...-z^" , for a = (ai, . . . ,a<i) e Z'' . 
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In the symbol notation, the subdivision step in (|l.ip is described by the identity 

{S,d){z)^d{z')a{z) , z'^ = {zf,zl...,zj). (1.2) 

The first factor on the right-hand side of (|1.2p refers to an upsampled version 
of the data d. 

Equation (jl.2[) can also be written using convolution with the so-called sub- 
masks of a. Let 

E={0,1}'' (1.3) 

be the set of representatives of Z'^/2Z'*, given by the vertices of the unit cube 
[0, 1]^, containing 

0== (0,0,...,0) and 1 = (1, 1, . . . , 1) . 

Then, the 2"^ submasks ag and their symbols ae{z) are defined by 

ae = (ae+2a)aez<i and a^iz) = ^ ae+2a , eeE. 

The standard decomposition 

a{z) ^^z"" a^{z^) 

eGE 

yields the equivalent form of the identity (|1.2p 

(5ad)(2)-^2^d(^2)a,(22). 

eeE 

This shows that a subdivision step is the result of convolving the input data d 
with each submask ae, which is a process of low pass filtering, followed by an 
interleaving process, i.e., upsampling and multiplication by z^, to produce the 
output data S^d. 

We say that the subdivision scheme Sa is convergent, if for any starting 
sequence d e £oo{'^'^), there exists a uniformly continuous function /d such that 

lim sup |(5:d)„-/d(2-'^a)| = 0, 

and /d 7^ for some initial data d. This is the notion of C-convergence, also 
referred to as uniform convergence, see where Lp-convergence for 1 < p < oo 
is also discussed. The necessary condition for these types of convergence is now 
known to be the so-called sum rule of order 1 referring to the submasks 

ae(l) = ^e+2a = 1 , eeE, (1.4) 
see [1, Proposition 2.1] and [2^, Theorem 3.1]. 
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Another notion of convergence is used in the hterature on multiresolution 
methods devoted to wavelet and frame constructions, for details see 0, Chap- 
ter 13], [13, Chapter 7], or [s^, Chapter 6]. There, the convergence is charac- 
terized by the properties of the infinite product 

na#(^^^), a#(^) = a(^), 

3=0 

if z is restricted to the d-dimensional torus. To state some of the properties we 
switch to the real variables ^ = ("^i, • • • Ad) via the transformation zj = e^™^\ 
j = 1, . . . , d. The set (|1.3p then transforms into the set 

Z = ZE = {e = e-"" : e e E} {-1 , +1}'^ (1.5) 

of the vertices of the cube [—1, +1]^*, and the necessary condition (|1.4p takes the 
equivalent form 

a(l) ^ 2'' and a(£) = for £ e Z' = Z \ {1} . (1.6) 



For this reason we call Z' the zero set, and the conditions in (|1.6p the zero 
condition of order one (Condition Z\). In the literature, both the conditions in 
(|1.4p and their equivalent form in (|1.6p are called the sum rules of order one. 

More generally, we also use the higher order sum rules following the nota- 
tion introduced and discussed in the survey paper 25[, see also the references 
therein: The mask symbol &(£) is said to satisfy the zero condition of order k 
(Condition Zf.)^ if 

a(l) = 2^^ and (i:'ja)(£) = for £ e Z' = Z \ {1} and |j| < fc . (1.7) 



2. Zero condition and the associated ideal 

In this section we show that condition Zi fully determines the structure of 
the symbol a.{z) of any convergent scalar subdivision scheme, see Theorem A. 
The statement of this theorem is therefore the first step toward the desired 
multivariate generalization of the Factorization Property of univariate schemes. 
Theorem A is also crucial for understanding the structure of the symbol &{z) 
satisfying condition Z^ for A; > 1. It also shows that the Factorization Property 
is replaced in the multivariate case by the decomposition (|2.2p . 

The generators q& in (j2.2p are the symbols of subdivision schemes whose 
limit functions le are the characteristic functions of the parallelepiped spanned 
by the column vectors of certain matrices 0. These functions are box splines of 
degree zero. In order to control the size of the mask, it is, thus, of importance 
to choose these box spline symbols appropriately, and to work in (|2.2p with as 
few summands as possible. Theorem 12 .61 shows how to meet these requirements. 

A modification of the generators q& is given in Theorem l2 . 1 1 1 and Theorem A, 
at the end of this section. The modified generators 5© have a useful algebraic 
property: they factor into d linear polynomials. 
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Since all convergent subdivision schemes satisfy condition Zi, we start with 
a characterization of the polynomial ideal 



I={peU'^ : p{e) = for e G Z'} , (2.1) 

and later, in Section 3, of its powers 

= {p e n"* : (£>>)(£) = for £ £ Z' , |j| < k] 

for k > 1. The main result of this section, whose proof is a consequence of 
Theorem I2.6[ states the following: 



Theorem A. The mask symbol of any convergent d-variate subdivision scheme 
iSa can be written in the form 

a(2;) = ^A0a0(2;)2'^g0(^) . (2.2) 


The sum runs over all unimodular dxd-submatrices & of (X^^^^ iX^^"^^ from 
(|2.6p . The polynomials q&{z) are defined in (j2.7p . a@{z) are Laurent polyno- 
mials satisfying cr©(l) — 1 , and A© are real numbers subject to A© = 1 . 

Remark 2.1. Since the 2'^q&{z) are the mask symbols of certain box splines 
of degree zero, the mask a is an affine combination of masks each of which 
originates from such a box spline convolved with some (smoothing) factor. Thus, 
Theorem A explains why such affine combinations were successfully studied 
before, see e.g., 

UMm for 

examples of bivariate and univariate schemes. 

From the point of view of algebraic geometry. Theorem A tells us that the 
system of box spline symbols g© in the representation ^2.2\i generates the ideal 
I. This connects our work to the papers 3^ 31 1 where the author has studied 
these ideals in detail and has characterized the ideal I in [sol. Proposition 4.1] 
using the generators 

I^<zl-l,zj-l,...,zj~l, (zi + 1)(Z2 + 1) . . . (zd + 1) > • (2.3) 
The characterization (|2.3p is a special case of his more general results, see also 



28l . Example 4] . It also follows from the fact that the polynomial ideal 

J ■.^{pen'^ : pie) = for £ e Z} 
of polynomials vanishing on Z is generated by the polynomials 

z\ ~ 1, ^2 — 1, . . . , Z^ — 1 . 

This result is stated in 0, Lemma 2.3] with an elementary and constructive 
proof. Consequently, in order to determine a set of generators for the quotient 
ideal I, it suffices to add to the generators of J the polynomial 

(Zl + l)(z2 + l)...(zd + l) 
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Z2 




Figure 1: The ideal X from the algebraic geometer's point of view. 



that vanishes on Z' , but does not vanish at 1 . 

In the two-dimensional case, guided by ideas from algebraic geometry, it 
seems natural to use the straight lines through the three points of Z', as shown 
in Figure [U in order to find other generators for the ideal I . This yields the 
following result. 



Proposition 2.2. For d — 2 , we have 

2" _ , 1 + 21 1 + Z2 1 + zi zi + Z2 1 + Z2 z\ + Z2 ^2 ^^ 

~ 2 2 ' 2 2 ' 2 2 ■ ^ ' ' 

Proof. By we have for rf = 2 

I^<zl-l, z| - 1, (21 + 1){Z2 + 1) > . 

Let us denote 

~ _ 1 + Zl 1 + 22 1 + Zl Zl + Z2 1 + Z2 Zl + Z2 

^ 2 2 ' 2 2 ' 2 2 ' 

then it suffices to show that the generators for I are contained in I and vice 
versa. 

To this end, note first that 

2 1 . 1 + Zl 1 + Z2 , . 1 + Zl Zl + Z2 

Zl — I = —4 h 4 ■ 



2 2 

1 + Zl 1 + Z2 , 1 + ^2 21 + 22 



1 = -4^— ^ — ^ + 4 

2 2 



/ I 1 \ / I 1 A - 1 + 21 1 + 22 
(Zl + 1)(Z2 + 1) =4— — 



8 



and conversely that 



= - (Zl + 1)(Z2 + 1) 





2 




2 




1 


+ zi 


Zl 


+ 


22 




2 




2 




1 


+ Z2 


zi 


+ 






2 




2 





= -(Zf-1) + -(Z1 + I)(z2 + 1) 
= i(^2-l) + i(^l + l)(^2 + l). 



This shows that X — I as claimed. 



□ 



Remark 2.3. The three functions in ()2.4p are remarkably close to the symbols 
of the three box splines of degree zero on the three-directional grid. Indeed, if 
we reflect the standard three-directional grid about one of the coordinate axes, 
i. e., use the grid spanned by the three vectors (J) , Cj*), and (li) , then the 
box splines of degree zero have the mask symbols 



1 + 21 1+22 



1 + 21 1 + 21/22 



1 + Z2 1 + 21/22 



Since we can write zi -f Z2 = Z2 (1 -I- Z1/Z2) , these are just the functions in (|2.4p 
up to an appropriate normalization. 

It follows that the result of Proposition 12.21 can as well use the usual box 
spline symbols: The reflection of the three-directional grid about the zi-axis, 
as mentioned before, corresponds to the variable transformation (zi , Z2) 
(zi , I/Z2) , and then (12. 4p becomes — again after a proper normalization — 



1 + 21 1 + 22 1 + 21 1 + 21 22 1 + 22 1 + 21 22 



(2.5) 



From the algebraic geometer's point of view, this amounts to replacing the 
straight line zi + Z2 = in Figure [T] by the hyperbola 1 + zi Z2 = . 



The identity (|2.5p is the statement of Theorem 12.61 for the case d = 2 . To 
be able to state this result in general, we first need to provide some additional 
notation. As before, we denote by E = {0,1}'^ the set of vertices of the d- 
dimensional hypercube, and let E' = E \ {0} . Collecting the elements of E' in 



the matrix Xn 



Xn 



/I 





Vo 







yields 

1 1 

1 

1 ■ 

: 

: 



1 

1 

: 1 



1 : 

y(2) 



1\ 



1/ 



(2.6) 



(k) 

where each column of the submatrix X^^ contains exactly k entries equal to 1 . 
Thus x'^P consists of the standard unit vectors Cfc , k — 1, . . . , d, while x'"?^ 
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contains the vectors ej + e^, j ^ k, etc. We treat the columns of Xd as 
directional vectors from which we build the box splines of degree zero, i. e., the 
characteristic functions of certain parallelepipeds. This means that we take any 
d columns from Xd to produce a square submatrix of Xd ■ With each such 
, we associate the normalized polynomial 



'7eW=nH^' (2.7) 



l + z' 

where 9 runs through the columns of , with q&{l) = 1 . 
Proposition 2.4. For any dxd-submatrix of Xd in (|2.6I) . we have 

q@ el det = 1 (mod 2) . 

Proof. By definition, G I if and only if 

VeeZ': gQ(£)= i±5!=o. 

ee& 

Due to = e^'^'"'^^^ , e £ Z', this is equivalent to 

Ve e E' : 30 e : e^6> = 1 (mod 2) . 

/. e., the map 
has the property 

Ve e E' : L{e'^) ^ 0^ (mod 2) . (2.8) 

As we are only interested in the parity of det , we may employ the ring 
homomorphism Z — >■ Z2 = Z/2 Z mapping even and odd numbers to their coset 
representatives and 1 , respectively. It extends naturally to a matrix ring 
homomorphism Z'^^'^ -> {'^2)'^^'^ and thus induces a linear map L : {^2)'^ —>■ 
(Z2)'' satisfying detL = det L . Since the set E is a complete set of coset 
representatives of the vector space (Z2)'* = Z''/2Z'' , property p.Sp is equivalent 
to _ 

ker(L) = {O"^} , 

i. e., L is a vector space automorphism and thus 

detL 7^0 inZa. 

In other words, 

det = det L ^ (mod 2) , 
which yields the claim. □ 

The next result, which is important for the proof of Theorem A, states 
that, conversely, the ideal I is generated by a subfamily of the elements q@. It 
is worth noting that the condition det = ±1 is equivalent to having an 
integer inverse and, thus, to the fact that its columns generate the integer grid 
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Lemma 2.5. For k = 1, . . . ,d , we have 
zl-l e < |g© : dxd-submatrix of [x'^J^ \ X 
Proof. For d = 1, the lemma simply claims that 



(1) i y{2)^ 



det© = ±1, efe e e| > . 



2 1 1+2; 
Z-1 G < > . 



which follows from z^ — 1 = 2(2; — 1) 



For d > 2 , it suffices to prove that 

Zk-le<{q@:&eU^/^}> for k = l,...,d, 



(k) 

where the family ' is defined to be 



= I© : dx ((i-l)-submatrix of (X^^ 



^^X^),dct 



;0iefc) =±l}. 



The proof is by induction on d . For d = 2 , we have (-X"^^^ i -X"^^^) = ( g ? i ) 
From this, we find for k = I 

"^"={(0.(0} 

and thus the lemma claims that 

-I ^ ^ 1 + 22 1 + Z1Z2 ^ 

Zl - I & < — - — , > 



which follows from 



^,_l = 2^ii±^-2 



2 

1 + 21 22 



2 2 

The claim for k = 2 follows by symmetry. 

For the induction step, we consider first k = d+1 and write 



(2.9) 



Zd+i - 1 = -Zd+i (zd - 1) ^ + {zd - 1) i±^^ + {zd+i - 1) ^ . (2.10) 
By the induction hypothesis, we can write 



Zd-'i-= ^ P@{zi,...,Zd)q@{zi,...,Zd) 



(2.11) 



(d) 



for certain polynomials ■ For e U^'''^ , we define 



©' 









( 














6 


, &" = 




6 




1 






1 


\o ... 


oy 




1,0... 


1/ 



and ' 



P0 6 

1 

0/ 
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where 





/I . 


0\ 


p = 





1 6 




oo 


■ 

■ 1/ 



It is easily seen that this yields 0' , 0" , 0"' £ Z^^+t^^ ^itli 



9©'(zi,...,Zd+i) =g©(zi,...,Zd 
q©"(zi, . . . , Zd+i) ==g©(zi,...,z<i 



+1 

1 + 



2 ' 

1 + Zd Zd+l 



and (7e'"(zi,...,Zd+i) = q&{zi, Zd-i, Zd+i) ^ ^^"^ ■ 
Using these identities, we obtain from p. lip 

-Zd+i{zd~l)^-^ ^ ^ -'Zd+ip&izi,...,Zd)q@'izi,...,Zd+i), 

1 + Zd+_ 

and, by replacing Zd by Zd+i in (|2.1ip 



{zd-l)'^-^-^Y^ ^ ^ P0(zi,...,Zd)g0"(zi,...,Zd+i), 



(Zd+l - 1) ^'-^ = X! P©(2^1> ■ • • )2d-l,Zd+l)(7©"'(2:i, • ■ • ,^:<i+i) • 



Substituting these three identities into (|2.10p yields a representation for z^+i — 1 
of the desired form. A corresponding representation for Zfe — 1 with k < d can 
be obtained by a cyclic permutation of the indices. 

This completes the induction. □ 

We would like to emphasize again that a very important consequence of the 
following theorem, for subdivision schemes, is stated as the main result of this 
section in Theorem A. 

Theorem 2.6. The ideal I is generated by the elements q@ , where are the 
unimodular dxd-submatrices of {X^^^ \ X^^"^) from (j2.6l) . 



Proof. One of the inclusions follows from Proposition l2.4l In the light of 
the other inclusion is an immediate consequence of Lemma 12.51 and the fact that 



(zi + l)(z2 + l)...(zd + l) = 2'^g0(2) for & = Id = x'^JK 

□ 

It is worth noting that the condition £ in Lemma 12.51 is equivalent to 
the divisibility of q& by l+z^ , which seems natural since l+z^ obviously divides 
1— z^ . Furthermore, the fact that each matrix in Theorem 12.61 contains a 
standard unit vector as its column is a consequence of the following result. 
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Lemma 2.7. Any dxd-submatrix & of satisfies 

det = (mod 2) . 

Proof. A dxd-submatrix of -X"^ contains in each column exactly two entries 
equal to one (and the others are equal to zero). For the linear operator L 
introduced in the proof of Proposition I2.4[ this implies 

L(l^) = 1^0 = 2 • 1^ = 0^ (mod 2) , 

1^ € ker(L) . From this, we may conclude that det L = in Z2 and, thus. 



I. e. 

det = det L = (mod 2) . 

Examples 2.8. We list the generators q& oil for low-dimensional cases. 

• For d = 1 , we have Xi = (1) , so the only submatrix is = Xi with 

q@(z) = ^- , 
and I is the principal ideal generated by this function. 

• For d = 2, we have 

X2 = 



□ 



which yields the following three generators for I 
'1 



01 = 

02 = 

03 = 



with 



with 



with 



^©1(2) = 
q&3 (z) = 



1 


+ Z1 


1 


+ Z2 




2 




2 ' 


1 


+ Z1 


1 


+ zi 22 




2 




2 


1 


+ Z2 


1 


+ 21 Z2 




2 




2 



as stated in ()2.5|) . These three functions generate X minimally in the sense 
that no two of them generate all of X . 

• In the case c? = 3, we have 











1 


1 







(: 


1 





1 





1 


li) 







1 





1 


1 





The submatrices with even determinants are obtained by selecting the 
columns (1, 2, 4), (1, 3, 5), (1, 6, 7), (2, 3, 6), (2, 5, 7), and (3, 4, 7), det = 0, 
and (4, 5, 6), det = —2 (compare Lemma [2771) . The remaining (3) ^ 7 = 28 
submatrices have determinant ±1 and thus describe elements of X . Restrict- 



(1) 



X^^^) as stated in Theorem yields (^) - 4 = 16 



ing ourselves to (X 
elements generating the ideal X . However, this system is highly redundant: 
using an algebraic manipulation program, wc found that it contains 64 sub- 
sets of 10 elements each which generate X minimally (in the sense that in each 
case, omitting any one of the 10 elements does not yield a set of generators 
anymore). 
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Geometric interpretation and modification. 

The columns of Xd describe the edges and the various diagonals of the d- 
hypercube. By Proposition l2.4[ the polynomial is an element of X if and only 
if the columns of span a parallelepiped with an odd d- volume. Theorem 12.61 
states that in order to generate I , it suffices to restrict ourselves to the edges 
efe and the 2-surface-diagonals + e^., j ^ k, of the hypercube only, and to 
consider only parallelepipeds with c?-volume equal to one. Lemma 12.71 shows 
that these parallelepipeds all have at least one unit vector as an edge. 

Proposition 12.21 leads to an alternative approach. For d = 2, the three 
vectors ei — (J), 62 — ("j*), and ei ~ e2 — { can also be seen as the edges 
of the 2-simplex. Since for d > 2 we only need to use the diagonals of the 2- 
dimensional surfaces of the d-hypercube, it seems natural to consider the edges 
of the d-simplex instead. 

To this end, we describe a different family of polynomials. For any dxd- 



submatrix of (X: 



(1) : v(2)^ 



define 



q&{z) = r0{z) with re{z) 



2 ' 

Zj + Zfc 



if = Cfc , 

ii 6 = Ej + Efc , 



Geometrically, this amounts for d = 2 to replacing the vector ( i[ ) by ( ) as 
third direction in the three-directional grid. Equivalently, this can be seen as a 
reflection of the grid about one of the coordinate axes (compare Remark 12.31) . 

Algebraically, this has the advantage that each factor rg is linear and there- 
fore that all the share the same total degree d . These functions still satisfy 
the analogues of Proposition l2.41 Lemma [2?5| and, consequently, of Theorem l2.6l 
and Theorem A. 



Proposition 2.9. For any dxd-submatrix of (x|^^'' \ x'"^^') 

q& el dct = 1 (mod 2) . 

Proof For e = (si, . . . , 6^) G Z' = {-1 , +1}'* \ {1} , we have 



(|2.6|) . we have 



- £k 



£k 



1 + ej/ek 



± 



2 2 2 ' 

and, therefore, q&{e) = ztq^^e) . Consequently, 



q@ el 



q&\ 







q&\ 



q@ ei. 



The claim follows by Proposition 12.41 
Lemma 2.10. For k — 1, . . . ,d , we have 



l-zl e < {q© : C [X'^p \ X^^) , det = ±1 , efe e 0} > 



(2)n 



□ 
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Proof. The proof is mutatis mutandis the same as that of Lemma 12.51 except 
that (|2.9p becomes 

^ 2 2 

and (|2.10p becomes 

^d+i - 1 = (zd - 1) ^ [Zd - 1) + (Zd+i - 1) ■ 

□ 

Theorem 2.11. T/ie ideall is generated by the elements q& , where are the 
unimodular dxd-submatrices of {X^^^ \ X^^^^ from p.6p . 

Proof. As m the proof of Theorem l2.61 using Lemma [2 . 1 01 instead of Lemma 1^751 

□ 



The interpretation of Theorem 12.111 as a property of subdivision schemes 
leads to the following modification of Theorem A: 

Theorem A. The mask symbol of any convergent d-variate subdivision scheme 
iSa can be written in the form 

a(2;) = ^ Ae cr© (2;) 2"^ ge (2;) , 

where a@{z) are Laurent polynomials satisfying a@(l) — 1 , and A© are real 
numbers subject to X)© = 1- ^"^^ runs over all unimodular dxd- 

submatrices & of (-X"!^^'' i^l^^'') from (|2.6p . 



3. Zero condition of higher orders and powers of I for d — 2 

In this section, we describe families of generators for the ideals I*^ , for A: > 1, 
and the implications on the representation of mask symbols from Theorem A. 
The case d = 2 is of special interest in subdivision. Therefore, most of this 
section is devoted to the study of this case, see in particular Theorem C below 
and the remarks at the end of Section 3.1. 

We start with a simple observation from ideal theory, namely, the fact 
that the product Ii • I2 of two ideals Ii — < : j — l,...,n}> and 
X2 = <{bk ■ k = l,...,m}> in a ring is generated by the pointwise prod- 
ucts of the corresponding generating sets 

Ii -12 ^ <{ajbk : j ^ 1,. ■■ ,n,k ^ 1,. . .,m}> . (3.1) 

Applying this to the construction of generators for I*^, we get the following 
immediate generalization of Theorem A. 
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Theorem B. A convergent d-variate subdivision scheme satisfies the condi- 
tion Zk if and only if its mask symbol can be written in the form 

j 

where Qj{z) are k-fold products of Laurent polynomials q& with unimodular 
dxd-submatrices & of Xd from (j2.6l) . <Jj{z) are Laurent polynomials normalized 
by <7j{l) — 1 , and Xj are real numbers subject to Aj = 1 . 

Note that the normahzations assumed in Theorem B imply a(l) = 2"^ . 

In the bivariate case, we can be much more specific and show that the 
generators for I*^ are the mask symbols of certain three-directional box splines. 
The latter have been studied thoroughly in 0, H, [23] , and their mask symbols 
have the form Ba^p^-y = 4:Bf ^ ^ , where 

«,A7eNo, (3.2) 

are the normalized box spline symbols satisfying ^^^^(1,1) = 1- The three 
indices a, /3,7 correspond to the multiplicities of the three vectors ei, 62 and 
61+62. In this notation, (|2.5p becomes 

^=<Sfi,0:Sfo4:<M>- (3-3) 

Note that the family of these (normalized) box spline symbols is partially or- 
dered and closed under multiplication, since 

Furthermore, they satisfy the following relation that we need later on. 

Lemma 3.1. For any given triple {a,l3,j) G Nq, the ideal generated by the 
three symbols B'^^-^ ^ ^ , B"^ ^ , and i?^^ ^_|_i is the principal ideal generated 

Proof. Since each of the symbols -B^^^ ^ ^ , B'^ ^ , and B*^ is a multiple 
of B'^ p ^ , we only have to show that the latter can be generated from the former 



3* 

three. To this end, we make use of the identity 



i(l-z2)Bfo,o(^i,^2) + i(l-^i)<i,o(^i,^2) + <o,i(^i'^2) = 1. (3.4) 



Multiplying both sides by -Bjf ^ ^ proves the lemma. □ 



In the light of p.3p . it is natural to expect that the higher powers of the 
ideal X are generated by box spline symbols of higher order, as the following 
result shows. 
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Theorem 3.2. In the bivariate case, the k-th power , fc G N, of the ideal I 
is generated by the set of three- directional box spline symbols 



{ ^t,l3,a' -^f,a,P^ -^a,l3,l3 



# 



0,1,. 



Proof. The proof is by induction on k . For k = I, the claim is just the iden- 
tity |33| . 

For the induction step, we write 1'^+^ = I ■ 2'' and apply (|3.ip . Using once 



more 
with 



and the induction hypothesis X — < Ife > yields I' 



k+l 



B 



# 



/3+l,/3 + l,a ' 
R* 

-°/3+l,;3,a+l ' 

D# 



# 



# 



;3 + l.Q+l,/3 ' ^q;+1.^+1,/3 

, B* 



/3+l,a,/3+l : 



a+l,/3,;3+l' 
B* 



a — 0,1, 

So we need to show that < I'^,^-^ > = < Ifc+i > where 



I3^k-a 



T i 0# p# R7 

ife+1 — ^5,(5,7' -°(5,7,(5' ^7,(5,i5 



7 = 0,1,, 



(5 = fc + 1-7 



To this end, it suffices to show that both I^+i C < > and C < Ik+i > ■ 
Firstly, note that the elements of Ife+i are the diagonal elements in the list 
-^'fe+i ^ith J — a and thus S = , with the only exception 7 = > 
This can only happen if A: = 2 ^ + 1 is odd, and then 7 = ^ + 1. For this value 
of 7 , the list Ifc+i contains only one element, viz., Bf_^^ . By Lemma [XTl 

this is generated by the three elements above the main diagonal in with 
k = 2e + l anda = ^, /3 = £+1 . 

Conversely, the diagonal elements of are all listed in I^+i . The indices 
of any non-diagonal element form a permutation of the triple {/S, (3+1, a+1) . 
The associated is a multiple of B'^ ^ which appears together 



f3 , except for the case 



with all index permutations in Ife+i for 7 = a+1 , S 
k-\-l 

. This can only happen \i k — 21 is even, and then a — t 

and {fi, (3+1, a+1) = {£,£+l,£+l) . But the associated elements appear in I^+i 
for 7 = ^ , (5 = f +1 . □ 



As in Section [21 Theorem 13.21 has an immediate consequence. Theorem C, 
for bivariate mask symbols satisfying the higher order zero conditions. This 
result is of great importance for studying the properties of existing subdivision 
schemes and also as a starting point for the construction of new schemes. 

Theorem C. A convergent bivariate subdivision scheme iSa satisfies the condi- 
tion Zj. if and only if its mask symbol can be written in the form 



a{z) 



E 



A 



(3.5) 



where ^ A^.^^-y = 1 , and the <ya,i3,-yiz) are Laurent polynomials normalized by 
the condition crct^^.^(l) = 1 . 
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Examples 3.3. We illustrate the result of Theorem [321 by explicitly listing the 
generators for small values of k . 

T r D# 0# D# 0# "1 

^2 — t -^2,2,0 ' -°2,0,2 I ^0,2,2 ' -°1,1,1 J ' 

T r D# R# R# R# D# D# 1 

^3 — t ^3,3,0 ' -°3,0,3 ' -°0,3,3 ' -°2,2,1 ' -°2,1,2 ' ^1,2,2 J i 

T r R# R# R# R# R# R# R# 1 

M — t -°4,4,0 I -°4,0,4 ' -°0,4,4 J -°3,3,1 ' ^3,1,3 ' "1,3,3 ' -°2,2,2 / • 



3.1. Further properties. 

To be able to show some further properties of the box spline symbols, im- 
plying the corresponding properties of the associated subdivision schemes, we 
need the following auxiliary result. As usual, we write 

an + m 

for mixed partial differential operators. 

Lemma 3.4. The partial derivatives of the box spline symbol ^ are given 
by 



=0 ^1=0 

m\ ^ l3+£\ f l+Z2\^+'^-''"'-lh-y-e\ / Z2~iy-'^-l 



X 



2m \m-jJ \ 2 J \ j 

j=0 



Proof. The identity 

1 + 21 22 _ 1 +21 1 +22 .21 - 1 22 -1 

2 ~ 2 2 2 2 

yields 

f=0 
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From this we obtain with the Leibniz formula 



«=0 

7 



«=0 ^j:=o 



EC) EC) 



Ec; 



{/3+£)! (l+22)'3+''-(™--') (7-^)! (22-1)^" 



(/3+^-(m~j))! 2"+-^ (7-^~i)! 2^-*^ 

^ E G) (e iiO (^) r;') 



Em!_ / 1+Z2 \ 
2»"Vm.-7/V 2 / 



2 / 

as claimed. □ 



Together with Theorem 13.21 this allows us to determine the maximal order 
of sum rules satisfied by a three-directional box spline symbol. 

Proposition 3.5. For any triple (a,/?, 7) £ Nq , the maximal k such that 
^f,/3,7 ^ ^'^ ^■^ d^ven by 

k = a + /3 + ^~ max{a, /3, 7} . 

Proof. Let (a, /?, 7) = 7r(a', /?', 7') where tt is a permutation such that a' < (3' < 
7' . It follows from Theorem 13.21 that 

^a,l3,j — -°7r(a',0',7') ~ ^7r(0,0,7'-/3') -°7r(Q',/3',^') ^ ^ 

for = a' + /3' = a + /3 + 7 — maxja, /3, 7} . 
On the other hand, applying p. 61) yields 

and {D^-''^^B*^^^)il,^l) =(-l)^^^0. 

This shows that B*u^(^X^+^ for 

fc = min{a+/? , a+7 , /3+7} = a + /3 + 7 — maxja, /3, 7} , 
and this completes the proof. □ 
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The sets of generators described in Theorem l3.2l are minimal generating sets, 
as the foUowing shows. 

Proposition 3.6. The set 1^ generating l'' is minimal in the sense that for any 
p ^ e Ifc , the reduced set 1^ \ {B"^ ^ ^} does no longer generate the ideal . 



Proof. We begin with an element of the form B* ^ ^ G Ifc . We find that 

(^^"'"^</5,,)(-l:-l) = ^^0, (3.7) 

and we claim that for all other elements Bt G Ik Ifc \ {B^ o a \ , we have 

5# V_l^_l)=0 for all (0,0) < (n,m) < (a,/3). (3.8) 

To this end, assume first that a < /3 , then we have 

[dMb*^J){-1,-1) = Q and ) J (-f , -f ) ^ , 

since in p.6p . we have (5 + 1— {n — i)>j5 — a>{) and therefore we also have 
I J = for zi = — f . For arbitrary a < j3 , consider (a',/?') ^ 

In case a' < a and thus f3' > /3 , we have 

(-1,-1)^0 (3.9) 

since in (13. 6p . we have /3' + £ — {m — j) > /?' — /3 > and therefore we get 
( — - — 1 = for Z2 = — 1 ; and also 

(i?("'™)i?|,„,^,)(-l,-l)=0 and (l?("^'")i?f,^,„,)(-l,-l) = 

(3.10) 

since in (|3.6|) . we have /?' + — (n — «) > /?' — a > and therefore we also get 

(— j =Oforzi = -l. 

In case a < a' < < /3 , p.9p holds since in p. 61) . we have a' + f— (n — i) > 

a' — a > and therefore ( — ^ — j = for zi = — 1 ; and also p.lOp 

holds since in p.6p . we have (3' + £ — (n ~ i) > jS' — a > and therefore 

f-^j = for zi = -1. So dSU is shown. 

But this implies that 

Bt0,p ^ <Ifc>, 

since otherwise, we could write 
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which, employing the Leibniz formula, yields 

= E EC) E0(^'-'^-'"^^.,.-. 

D# ^ n— m— 

and this contradicts p.7p and p.Sp . 

By symmetry in zi and Z2 , it follows that also 

Bto.,pi<lk\{B*^^p}>. 

It remains to show that 

This can be achieved by employing a directional derivative and considering 
mixed derivatives of the form 

KdzJ \dz2J \dzi dz^J 

Along the same lines as above, one shows that 

(i?(-°^«<,J(i,-i) = (-iri^^o, 

but that for all other elements ^~ G Ifc \ {^f p , 

^£)(»,o/) ^^(1^-1)^0 for aU (0 , 0) < (n, £) < (a, /3) . 
Alternatively, we may employ the coordinate transformation 
{ ^ (^^ '^'^\ .i=> f^^\ ^ f "^^^ 

which yields 

and thus allows us to reduce this case to the first one also. □ 

The characterization of convergent bivariate subdivision schemes established 
in Theorem C opens a way for their systematic study. We only point out the 
following facts: 
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• The set 1^ of generators for I is symmetric in the sense that it is invari- 
ant under an interchange of the two variables, and that the indices of the 
generators can be permuted arbitrarily. 

• For even fc, the fc/2-th power of (the symbol of) the Courant hat function 
^1,1,1 appears in I/j . 

• Most interesting for us, however, is the fact that the smoothness of these gen- 
erators matches perfectly with the order of polynomial reproduction. More 
generally, the box spline symbol -B^^ ^ is an element of 

L^^^i) C C^'^"^) for K = a + ^ + 7-max{a,^,7} = fc, 

see This means that for each three-directional box spline, smoothness 
and polynomial reproduction match in the same way as in the univariate 
case. This property is not necessarily preserved under taking combinations 
as described in Theorem C, however, as we shall see below for the four- 
directional box splines. 

• It may also be of interest that the total degree of the box spline with the 
symbol B*^ij^^ is 

Ai = a + /3 + 7-2, 
which for the generators B'p ^ ^ , B'p ^ ^ , and B'^ ^ ^ with a < P yields 

^ = a + 2^-2 = 2fc-a-2 = 2K-a-2. 

3.2. Modification. 

As in Section [21 we can replace the q& by the q& , and all results still hold. 
More precisely, this means replacing B"^ ^ ^ by 

For example. Lemma 13.11 holds as stated, and in its proof, we only have to 
replace by 

■^fo,o(^l'22)+i?(f^i^o(^l'^2)-S(f^o,l(^l>2;2) = 1- 

Also, Theorem 13.21 together with its proof. Examples 13.31 and Theorem C are 
still valid. Lemma now runs as follows. 

Lemma 3.7. The partial derivatives of the modified box spline symbol ^ 
are given by 

-i:o(i:Scx^)°""""'0(l)"*) (3.") 
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Proof. Writing 

1=0 

yields by the Leibniz formula 



2-^ \U dz{' \ \ 2 ) \2 ) j dz^ \\ 2 ) V 2 / 



\i) {a-(n-i))\ 2° H-i)\ ' 



EG) EC 



X 



i/ {a-(n-i)y. 2° {e-iy. 2^ 

/3! (l+22)'^~('""^' (7-^)! zl"'^"-'' 



EC 



_ (/3-(m-j))! 2" (7-^-i)! 2^-' 

^EO(E^C)(ii^)'"""""{:)(?)' 



X 



^ 2™ Vm-1 
J=0 



as claimed. □ 



This implies that Proposition 13.51 remains valid, and in its proof, only the 

values of the listed nonzero mixed partial derivatives need to be multiplied by 

(— 1)^ . Also, Proposition 13.61 still holds as stated: in its proof, once more the 

signs of the nonzero derivative values change, and in the last part, we need to 

do 

consider the directional derivative — instead of h t: — ; alternatively, 

ozi 022 azi az2 

we can employ the coordinate transformation 

wi\ ^ /'zi/z2\ fzA ^ /'Wi/W2 

W2J \1/Z2) \Z2j \1/W2 

which yields 



Another approach is the use of the coordinate transformation 
W2J U2/ \}/w2 
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which was already mentioned at the beginning of Section [51 This yields 

which shows that all the results for the ^ ^ hold for the -B^^ ^ also. 

3.3. Four- directional box splines 

The refinement mask symbol of a four-directional box spline has the form 

This uses both alternatives, 61 + 62 and 61 —62, for the third direction discussed 
above. We can rewrite the above as 



Blp,,M^z^) = i^[} + iJ Bl^Jz,,Z2), (3.12) 

illustrating the well-known fact that any four-directional box spline is indeed 
a special convex combination of the shifts of some three-directional box spline. 
This convex combination uses the normalized binomial weights Xs^e — 2"'' (^) , 
i = 0, . . . , 6 . The representation in p.l2p is not optimal, however, if one tries to 
determine the maximal k such that B'j^ ^ ^ g G . Instead, we use the identity 
zi + Z2 = {1 + zi) {1 + Z2) — (1 + zi Z2) which in terms of the normalized box 
spline symbols reads as 

<0,0.l(^l>^2) = ^ ^ = ^ (2B*,^,{Z,,Z2)~B*,^,{Z,,Z2) 



This yields 



(zi,Z2) = <^,^(^i,^2) ^ (2i?# ^o(zi,Z2)-<o,i(^i'^2))' (3-13) 



B* 



s 



f=0 



Proposition 3.8. For any quadruple (a, (3, 7, 5) £ Ng , the maximal k such that 
Bt,i3,-y,5 S ^r'' is given by 



fc = a -I- /3 + 7 + (5 - max{a , /3 , 7-|-(5} . (3.14) 
a,/3,7,5 e if and only if B*^ ,^^^ 



other words, B'^ ^ ^ ^ G l'^ i/ flKC? on^j/ i/ B^ ^ G l'^ . 
Proof. By Proposition 13.51 we may conclude from p.l3p that B'^ „ g for 



min [a + l3 + -i + 5 + l~ max{Q;-f^ , P+t , -i+6-t}] 

£— 0,...,(5 

^ rnin ^{min{/3+7+(5 , a+7+(5 , a+(3+2i}) 

niin{;3+7+(5 , q:+7+(5 , a-l-/?} 

a + P + ^ + S — max{a , j3 , ^+5} . 
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Thus it remains to show that B^^ ^ s ^ X^^^ for this k . Equivalently, we show 
that 

~Bt,.M^^^) - ^/ - i}^T (^)' (^)' (^)' 

satisfies i?^^ ^ ^ ^ 1*"+^ for this k . 

To this end, we proceed along the hnes of the proofs of Lemmata 13.41 and 13 . 71 
and Proposition [3?5] Writing 



=EOE0m°"(^)^"(l)' 

f=0 j=0 




and, therefore, 

and (1 , -1) = (-l)-fteg^ ^ . 

This shows that B*^ ,^ g ^ X'^+i for 

k = min{Q!+/3 , a+j+6 , /?+7+(5} = q; + /3 + 7 + (5 — maxja , /? , 7+(5} , 
and this completes the proof. □ 

Remark. The proof of Proposition 13.81 establishes the connection between the 
order k in (|3.14l) and the smoothness of the (a, /3, 7, (5)-box spline, which can 
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be determined using the results in The function is an element of C 
C(«-2) for 

K ~ min{Q!+/3+7 , a+/3+5 , , /3+7+5} = q;+/3+7+(5— max{a , /? , 7 , (5} , 

while the order k (degree fc — 1) of polynomial reproduction is given by 

k — miii{a+/3 , a+j+6 , {3+^+6} = a + /3 + j + S — maxja , /? , j+S} . 

Corollary 3.9. For the four- directional {a, P,^,5)-hox spline, we have k> k , 
and 

K > k 7 + (5 > max{a, /?} and min{7, (5} > . 

Proof. By the above remark, 

K — k = max{a , /3 , 7+5} — max{a , /3 , 7 , 5} > . 

In the case max{a, /?, 7+5} — max{a, /?} , we have max{7, 6} < j+S < niaxja, /?} , 
i. e., also max{a, /3, 7, 6} — max{a, /?} and therefore k = k . Otherwise, we have 
7+(5 > max{ck,/3} , and then 

K — k — 7+5 — maxja, /?, 7, S} = min{7+5— a , 7+5—/? ,7,5} 

where the first two elements are positive, so in this case, 

K — k — <^=^ min{7, 5} = . 

All in all, this yields that 

K — k — 7+5 < max{a, /?} or min{7, 5} = , 

which proves the claim. □ 

The condition min{7, 5} = is worth a closer look. For 5 = , we are in the 
standard thrcc-directional case, while for 7 = 0, we have 

-Sf,/3,0,5(^l'^2) = B*p^s{zi,l/Z2) = z:^^ B*pj,{zi,Z2) , 

so these are three-directional splines on the refiected grid. 

Examples 3.10. We list a few standard examples together with the decompo- 
sitions of their mask symbols according to p.l3p and to p. 51) . 

The (1, 1, 1, l)-box spline, known as the Zwart-Powell element, has the mask 
symbol 

4-Bfi^i_i(zi,Z2) = ^ (2 B#2,i(^i' ^2) -5^1,2(^1,^2)) 
= 1 5 4B{^i i(zi,Z2) . 
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We find fc = 4 — 2 = 2, so the associated subdivision scheme reproduces poly- 
nomials of total degree up to one; but k = 4— l = 3,so the function is in 
Lg^ C C(i) . 

The (2, 2, 1, l)-box spline has the mask symbol 

= 2 • 1 • 4 <3.i(zi, Z2) + (-1) • f ■ 4i?#2,2(zi, Z2) . 

Here, we have = 6 — 2 = 4 , telling us that polynomials of degree up to 3 
are reproduced, and also k = 6 — 2 = 4,j. e., the function is an element of 
Lg^ C C(2) . 

More interesting are the higher order four-directional splines. For example, 
the (4,4, 1, l)-box spline has order of polynomial reproduction fc = 10 — 4 = 6 
and K = 10 — 4 = 6 . Its mask symbol can be represented as 

4<4,l,l(^l'^2) - ^ (255*5,1 - <4,2(^1>^2)) 

-2.1.4 B5*5,l(^l> ^2) + (-1) • f ■ 4i?#4,2(zi, Z2) . 
^2 22 

Remark 3.11. The vector case, i. e., the case when the matrix mask A e £q^'^{Z'^) 
is a finitely supported sequence of s x s-matrices indexed by Z**, is more intricate. 
The formulation of the zero conditions (jl.7p for the mask symbol of multivariate 
vector subdivision schemes depends greatly on the so-called rank of the scheme, 
see [fi, .12, 26]. Such a formulation, see, e. q., [23, does not allow us to read off 
the properties of the entries of the matrix Laurent polynomial A{z) directly. It 
is possible, though, to use a slight modification of the transformation in [22| to 
obtain the matrix sequences T, T"^ € t^'^'iZ'^) such that 

A{z) ^ T"''{z'^) ■ A{z) -Tiz) 

satisfies the zero conditions of a form that makes the structure of some of the 
entries of A(z) more evident. However, this is a topic for further investigations. 

4. Examples 

In this section we illustrate the result of Theorem C with some examples. We 
would like to emphasize that this result does not only simplify the study of the 
properties of subdivision schemes, but also yields a way for enhancing certain 
properties of existing schemes by combining them appropriately. 

In the following, the set 1^ of the box spline symbols is the set of generators 
for l'^ as in Theorem 13.21 

In the masks displayed below the boldface entry at bottom-left position 
refers to the index (0, 0). This assumption is not really important, but as stated 
above already, we prefer to shift masks so that they are supported in the positive 
quadrant and have polynomial symbols. 
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4.1. A bivariate interpolatory scheme 

Interpolatory schemes are characterized by the fact that one of the sub- 
masks is a (5 sequence, or equivalently, one of the subsymbols is identicahy one. 
Theorem C allows us to present a systematic way for creating interpolatory 
schemes from our lists of generators by equating the coefficients of their affine 
combinations and normalizing them appropriately. 

To provide just one such example, consider the interpolatory scheme studied 
in 23|, Example 2], a bivariate version of the univariate four-point interpolation 



scheme given in [18|. Its mask is 



1 

32 








-1 


-2 


-1 




























-1 





10 


18 


10 





-1 


-2 





18 


32 


18 





-2 


-1 





10 


18 


10 





-1 




























-1 


-2 


-1 





0/ 



The scheme reproduces polynomials up to degree fc — 1 = 3, whence a.{z) G X'', 
and a representation of a.{z) in terms of three-directional box splines from the 
list I4 is given by 



z^z^ a(zi,Z2) 



2* b|^4 o(zi, Z2) - 2 (z^ + Z2') Bl^^^{zi,Z2) + 2^ (1 + zi + Z2) B|^3^i(zi, Z2) 



2 
1 

2 I ,2 



2^ D# 



4 { -4 (^1 , ^2) - <2,2 (^1 , ^2 



3 3,3 



l(zi>22)} . 



From the second line, the weights A are recognized as —4, —1, and 6, and the 
normalized cr-symbols are 

zl + zl I + Z1+Z2 

1 , , and , 



respectively. 



4-2. The butterfly scheme 

The butterfly scheme has been studied in |21| and |23j, Example 5] . Its mask 
is given by 



/ 











-1 


-1 


o\ 








-1 





2 





-1 





-1 


2 


8 


8 


2 


-1 








8 


16 


8 








-1 


2 


8 


8 


2 


-1 





-1 





2 





-1 








^0 


-1 


-1 
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It is an interpolating scheme, and reproduces polynomials of degree A: — 1 = 3. 
The representation of the mask symbol in terms of three-directional box spline 
symbols from the list I4 is given by 



a(zi, Z2) 



^^^^^-S*3,i(^i'^2) - 2 zaSf 13(21,22) - 221^^3 3(21,22) 



13 



= 4|26- 
l + 21 + ^2 n# / x\ 
- 21 Bl2 2[Zl,Z2)j 

= ZIZ2B*.2 2{Z1,Z2) 

- 2 2lS]^3 3(21, 22) - 2 22^1^1 3(21,22) - 2 2122^1^3 1(21,22)} ■ 

We see that the generators are all multiples of Bf-^ 1(21, 22). This tells us that 
the symbol can be factorized as 

a(2i , 22) = 11 (21 , 22) b(2i , 22) , 

a fact noticed in [iTl . We would like to emphasize the following properties of the 
butterfly scheme. Firstly, a simple computation yields that the symbol h{z) does 
not define a convergent subdivision scheme, although each of the summands in 
h{z) by itself does correspond to a convergent scheme. Secondly, butterfly is 
an interpolatory subdivision scheme, but none of the summands in the affine 
combination above possess this property. 

^.3. A convergent scheme 

The symbols presented in the above examples all possess a property that is 
very important for their regularity analysis: they are multiples of one specific 
box spline symbol of type B'^ ^ ^ . The regularity analysis of such schemes is 
given in [3, Section 4.3]. The type of factorization used there, however, is a 
very special situation which does not generally hold for convergent schemes. 

A very simple example that comes to mind is the symbol given by 

a(2l,22) = 4 1^5^1,0(^1,^2) + ^5*1,2(21,^2)! = 4 C(2l,22) 

, . 1 1 + 21 1 /^l + 2l22 

With C(2l, 22) - + - 

By Theorem 13.21 the symbol &{z) is in I, but none of the generators from the 
list Ii divides the symbol. 

In order to check the convergence of this scheme, we study the properties 
of the so-called difference scheme 5b , see [3, [sij , with the matrix mask symbol 
B(z) satisfying 

i\~ l!)' - ^i^) ■ 
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One possible such B{z) is given by 



r,/, ^ N /^bii(zi,Z2) bi2(zi,Z2) 
^ ^'~[h2l{zi,Z2) b22(zi,Z2) 

with bii(zi,Z2) = i (zizf - + zizf + z| + 4z2 + 2) , 

hl2{zi,Z2) = , 

^2l(zi,Z2) = ^ (Zl2;2 - Zi - Z2 + 1) , 

and h22{zi,Z2) = - (^i^l + 2ziZ2 + 2zi + 3) . 

To check the convergence, we have to verify that the vector subdivision scheme 
5b converges to zero, see [1, The symbohc calculations yield ||5|||oo < 1 
and, thus, that 5a is C-convergent. 

It is worth noting that in this example the two building blocks, with sym- 
bols 4 Bf^ 0(21,2:2) and 4 Bf^ 2(^11^2), are not the symbols of C-convergent 
subdivision schemes, while the combination yields C-convergence. 

We also refer to the constructions in [lO|, where the convex combination 
of a four-directional, zero order box spline and a C^'^-quadratic box spline are 
used to obtain the so-called GP pseudo-quadratic box spline. This example 
shows enhancement with respect to linear independence of the translates, at 
the expense of reduced joint smoothness. 
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